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FOREWORD 


The  reliability  of  a  complex  system  frequently  enters  into 
important  decisions.  The  decision-making  process  may  be 
difficult  when  the  lifetime  distribution  of  the  particular  system 
is  known;  it  becomes  more  difficult  when  the  lifetime  distribu¬ 
tion  is  unknown.  Such  maybe  the  case  when  a  new  type  of  equip¬ 
ment  is  placed  in  operation  on  a  large  scale.  Provisioning  for 
replacements  should  be  based  on  the  lifetime  distribution  of  the 
equipment,  but  a  reliable  estimate  of  this  distribution  is  usually 
not  available  ur.til  a  great  deal  of  failure  information  has  been 
obtained.  In  such  a  situation  it  is  the  aim  of  the  Bayesian  ap¬ 
proach  to  enable  all  available  data  to  be  incorporated  into  the 
decision-making  process. 

This  paper  first  discusses  some  characteristics  that  are 
desirable  in  the  ckss  of  lifetime  distributions  to  be  examined 
by  Bayesian  methods  and  shows  that  the  Weibull  distribution 
possesses  some  of  these  characteristics.  The  remainder  of  the 
paper  treats  the  case  in  which  the  shape  parameter  of  the 
Weibull  distribution  is  known  and  presents  prior. posterior, and 
preposterior  Bayesian  analyses. 


Nicholas  M.  Smith 

Head,  advanced  Research  Department 
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Use  of  the  Weibull  Distribution 
in  Bayesian  Decision  Theory 


ABSTRACT 


The  Weibull  distribution  is  useful  in  analysing  the  probabilistic 
lifetimes  of  many  electrical  components  and  complex  systems.  It  is 
attractive  for  Bayesian  declsion-maldng  because  its  right-band  ci  nu- 
lative  function  is  of  an  exponential  form  which  allows  all  life-teat  data 
to  be  easily  incorporated  into  the  decision-making  process.  Un¬ 
fortunately  no  natural  conjugate  prior  distribution  exists  if  both  the 
shape  and  scale  parameters  of  the  Weibull  distribution  are  assumed 
to  be  unknown.  If  the  shape  parameter  is  assumed  known,  however, 

Bayesian  analysis  becomes  little  more  difficult  than  for  the  expo¬ 
nential  distribution,  a  special  case  of  the  Weibull.  Prior,  posterior, 
and  preposterior  analyses  are  given  for  the  case  of  known  shape 

parameter.  In  connection  with  preposterior  analysis  several  sampling  t 

plans  are  discussed.  The  paper  concludes  with  an  analysis  of  a 
problem  in  optimal  sampling. 


I 


1.  INTRODUCTION 


Many  decisions  Involve  the  reliability  of  electrical  components  and  com¬ 
plex  systems.  The  decision-making  process  is  often  difficult  when  the  lifetime 
distribution  of  the  particular  component  or  system  is  known;  it  becomes  more 
difficult  when  the  lifetime  distribution  is  unknown.  Such  is  the  case,  for  ex¬ 
ample,  when  a  new  type  of  equipment  is  placed  in  operation  on  a  large  scale. 
Provisioning  decisions  should  be  based  on  the  lifetime  distribution  of  the  equip¬ 
ment,  but  a  reliable  estimate  of  this  distribution  is  not  usually  available  until 
a  considerable  amount  of  failure  information  has  been  obtained.  This  paper 
first  briefly  discusses  the  benefits  that  may  be  obtained  by  applying  the  Bayesian 
approach  to  decision  problems  involving  the  reliability  of  systems.  Section  2 
then  discusses  some  desiderata  of  the  class  of  lifetime  distributions  to  be  ex¬ 
amined  by  Bayesian  methods  and  shows  that  the  Weibull  distribution  possesses 
some  of  these  desirable  characteristics.  Section  3  presents  a  Bayesian  analysis 
of  the  Weibull  distribution  when  the  shape  parameter  is  assumed  known.  This 
assumption  is  carried  through  the  remainder  of  the  paper.  Experimentation 
and  preposterior  analysis  are  discussed  in  Sec4,and  Sec  5  illustrates  the  ap¬ 
proach  with  an  optimal  sampling  problem. 

In  the  Bayesian  approach  we  assume  that  the  probability  density  function 
f(x)  of  the  useful  lifetime  7  (a  tilde  indicates  a  random  variable)  of  a  system 
is  a  member  of  a  class  of  density  functions  indexed  by  a  parameter  vector  0. 
Thus  we  write  f(x |0).  We  assume  that  the  value  of  0  which  applies  to  the  par¬ 
ticular  system  under  study  is  unknown  and  that  the  decision-maker  treats  0  as 
a  random  variable  and  expresses  his  Judgment  about  0  in  the  form  of  a  proba¬ 
bility  distribution  on  0.  We  write  the  probability  density  function  of  S’  as  f(0|ft), 
ft  indicating  the  experience  of  the  decision-maker.  By  using  the  density  func¬ 
tions  f(x|0)  and  f(0  |ft)  and  his  utilities  for  consequences  associated  with  future 
events  involving  x  and  S’,  the  decision-maker  can  choose  among  various  courses 
of  action.  If  some  evidence  z  (failure  data)  is  accumulated,  the  decision-maker 
revises  his  probability  distribution  on  <Tby  Bayes'  theorem: 

m Ms,  (1) 

/r<i|0.ft)f<0|ftW« 

Here  ilz  |0,  ft)  is  the  likelihood  of  the  evidence  z  conditional  on  0  and  ft . 

The  decision-maker  will  now  choose  among  various  courses  of  action  by 
employing  the  revised  or  posterior  density  function  f(8  |r,  ft).  Indeed,  whenever 
a  decision  must  be  made,  the  total  evidence  z  accumulated  up  until  that  time 
can  be  incorporated  into  the  decision-making  process  via  the  new  density  func¬ 
tion  f(0  |z,  ft).  This  is  one  of  the  main  benefits  of  the  Bayesian  approach. 
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The  notion  of  a  sequence  of  decisions  raises  the  possibility  of  deciding  to 
delay  a  decision  until  more  evidence  is  available.  But  Just  how  much  evidence 
should  be  collected  before  a  terminal  decision  is  made?  The  decision-maker 
can  decide  this  by  using  the  density  functions  f(x  |8)  and  f(0|£)  and  his  utilities 
lor  consequences  associated  with  future  events  involving  7,  %  and  the  evidence 
7  (which  is,  a  priori,  a  random  variable).  This  is  another  major  benefit  of  the 
Bayesian  approach  to  decision-making. 

The  mathematical  foundations  of  Bayesian  decision  theory  are  set  forth 
in  Raiffa  and  Schlaifer.1  For  further  discussion  of  the  Bayesian  approach  to 
decision  problems  involving  the  reliability  of  systems  the  reader  is  referred 
to  Briggs,8  Howard,*  and  Martel.'1 


2.  DESIDERATA  FOR  BAYESIAN  ANALYSIS  OF  RELIABILITY  PROBLEMS 

We  first  list  some  desirable  properties  of  the  density  functions  f(x|8). 

(1)  f(r|8)  must  be  the  density  function  of  a  nonnegative  random  variable. 

(2)  Depending  on  0,  f(x|0)  must  assume  a  fairly  wide  range  of  shapes. 

(3)  The  likelihood  -t  ( 2 1 e,  6)  should  have  a  relatively  simple  form. 

The  first  property  is  clearly  necessary.  The  second  is  desirable  if  the 
analysis  is  to  cover  a  wide  range  of  reliability  problems.  For  example,  the 
exponential  distribution  given  by  f(x  1 8)  =  e~0x  assumes  that  the  failure  rate 
does  not  change  with  time,  a  somewhat  restrictive  assumption. 

At  any  time  the  total  evidence  z  will  consist  of  the  information  that  r  sys¬ 
tems  have  failed  after  operating  for  times  Xj  ,  .  .  .  ,  xr  and  that  it  -  r  systems 
have  operated  for  times  xr+I  ,  .  .  .  ,  x„  without  failing.  We  assume  that  systems 
fail  independently.  If  we  define 


G(ijO)  >  r/WWy  . 

z 

the  likelihood  of  all  this  evidence  is 

f(*|0.6) .  n(  G(x,|t«J  . 

Unless  this  likelihood  has  a  relatively  simple  mathematical  form  the  compu¬ 
tations  required  by  Bayes’  theorem  become  very  complicated.  For  example, 
the  gamma- 1  distribution  with  density  function 

f(x|r.y)  .  y'z'-lt-yx/rir), 

where  =  (7,  y),  does  not  yield  a  closed-form  expression  for  G(x|r,  y)  unless 
r  is  an  integer,  and  even  then  G(x|r,  y)  is  a  sum  of  r  terms. 

A  density  function  which  has  the  three  properties  given  above  is  that  of 
the  Weibull  distribution  (see  Qureishi*  for  references): 

fw(xU,  oO  -  Aax0_,«~A,a. 


Gv(i|a.  a)  «  f  A,a. 


1 

f 

t 


d 

{ 
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The  parameters  a  and  A  are  shape  and  scale  parameters,  respectively.  The 
likelihood  l(z\6,&)  becomes  In  this  case 

fd|A.  a)  -  a'/Cr,  .  . .  t,)*'1  exp  ^  1“  J  .  (3) 

We  have  just  discussed  some  properties  that  are  desirable  in  the  density 
function  /(xj0).  There  are  also  some  properties  that  are  desirable  in  the  prior 
density  function  f(9  j  S ).  These  are  discussed  by  Raiffa  and  Schlaifer,1  and  we 
list  three  of  them  here. 

(4)  It  should  be  fairly  easy  to  ascertain  the  posterior  density  that  results 
from  a  given  prior  density  and  given  evidence  z. 

(5)  f(d  |  S)  should  be  a  member  of  a  closed  family  /*  of  density  functions 
so  that  the  posterior  density  f(0|z,  S)  is  also  a  member  of  f*. 

(6)  The  expectations  of  some  simple  utility  functions  with  respect  to  any 
member  of  f *  should  be  expressible  in  convenient  form. 

Unfortunately  it  does  not  seem  possible  to  find  a  family  f*  that  possesses 
these  properties  when  the  parameters  A  and  a  of  the  Weibull  distribution  are 
both  assumed  to  be  unknown.  This  is  basically  because  sufficient  statistics  of 
fixed  dimensionality  (see  Raiffa  and  Schlaifer1)  for  the  evidence  z  do  not  exist 
in  this  case.  Appendix  A  briefly  examines  a  joint  density  function  for  (A,  a) 
that  has  property  4  and  has  property  5  in  a  weak  sense  but  does  not  possess 
property  6. 

It  thus  appears  that  the  Weibull  distribution  with  both  parameters  unknown 
is  not  amenable  to  useful  Bayesian  analysis. 


3.  BAYESIAN  ANALYSIS  OF  THE  WEIBULL  DISTRIBUTION  WITH 
KNOWN  SHAPE  PARAMETER 

In  the  remainder  of  this  paper  we  shall  assume  that  f (x |0)  is  given  by  the 
Weibull  density  function  (expression  2)  and  that  the  valuejof  the  shape  param¬ 
eter  or  is  known.  In  this  case  only  the  scale  parameter  A  need  be  treated  as 
a  random  variable.  The  following  have  been  offered  as  justification  for  this 
assumption: 

(1)  For  some  items,  such  as  vacuum  I  »bes  (see  Kao*),  an  appropriate 
value  of  or  may  be  known  from  previous  test  evidence. 

(2)  The  same  assumption  is  made  in  standard  Weibull  life-testing  pro¬ 
cedures  (see  Ref  7). 

(3)  Previous  Bayesian  analysis  of  the  exponential  distribution  (see  Raiffa 
and  Schlaifer,1  Briggs,*  and  Martel*)  in  effect  makes  the  assumption  a  -  1.  It 
i«  no  less  valid  to  assume  a  value  of  or  other  than  1  in  an  appropriate  situation. 

Natural  Conjugate  Distribution 

We  need  to  find  a  family  f*  of  density  functions  that  possesses  properties 

4,  5,  xnd  6  of  Sec  2;  one  such  family  is  that  of  natural  conjugates  to  the  sample 
likelihood  (expression  3).  With  a  known,  a  kernel  of  the  likelihood,  i.e.,  a 
factor  which  depends  on  A  (see  Raiffa  and  Schlaifer,1  p  30), is 

A'<-Ay.  (4) 


5 


where 


(5) 


When  X  is  a  random  variable,  the  natural  conjugate  of  expression  4,  obtained 
by  treating  expression  4  as  a  kernel  of  a  density  function  for  X  (see  Raiffa  and 
Schlaifer,rp  47),  is  the  gamma- 1  density  function. 


fy,  y* 


v'x 

~m — 


0  <  A  <  -. 

0  <  r,  y  <r  -  . 


(6) 


Wilson*  han  previously  used  expression  6,  in  a  different  context,  as  a  prior 
density  function  when  the  scale  parameter  X  of  a  Weibull  distribution  is  unknown. 


Prior  and  Posterior  Analysis 

Suppose  the  decision- maker  must  choose  an  act  a  from  a  set  A  of  possible 
acts  and  his  terminal  utility  (terminal  in  the  sense  that  no  experiment  is  to  be 
performed)  for  an  act  a  and  particular  value  X  is  u(  (a,  X).  He  will  choose  an 
act  a'  such  that 


E^  u,  (a',  X)  *  max^  E^  ii,  (a.  A) . 

The  notation  F'  indicates  that  the  expectation  is  taken  with  respect  to  the 
gamma- 1  prior  distribution  on  X. 

If  a  gamma- 1  prior  distribution  with  parameters  r'  and  y'  is  assigned  to 
X  and  an  experiment  t  yields  an  outcome  z  (evidence  z )  with  statistics  r  and  y, 
the  posterior  distribution  on  X  is  gamma- 1  with  parameters  r  and  y",  where 

t"  •  r'  •  r,  y~  *  y'  ♦  y  .  (7) 

The  decision-maker  will  now  chooae  an  act  or  such  that 

Ej|,M((o",A)  •  maigE^i ,!«,(«,  A)  . 

Here  EJ(l  denotes  expectation  with  respect  to  the  posterior  distribution  on  X, 
which  is  conditional  on  z. 


4.  EXPERIMENTATION  AND  PREPOSTERIOR  ANALYSIS 

When  the  decision-maker  is  contemplating  a  specific  experiment  he  must 
consider  the  net  value  of  that  experiment.  He  will  generally  have  a  utility 
function  u(e,  z,  a,  X)  defined  for  each  combination  of  experiment  e,  outcome  z, 
subsequent  action  a,  and  particular  value  X.  Before  experiment  t  is  performed 
the  outcome  z  is  a  random  variable,  and  so  the  overall  (expected)  utility  of  ex¬ 
periment  e  is 

c  E,j(  m»*#  Ej[i  t  u(e,l.a.  A)  .  (8) 
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Here  Ei|<  indicates  expectation  with  respect  to  the  distribution  of  7.  When  a 
number  of  experiments  are  being  considered,  the  final  step  in  preposterior 
(before  an  experiment)  analysis  is  the  finding  of  an  optimal  experiment  e*.  i.e., 
one  such  that 

-  m»xf  u'it ) . 

It  is  clear  from  Eq  8  that  the  distribution  of  7  is  important  for  the  suc¬ 
cessful  completion  of  preposterior  analysis.  In  the  present  case  the  statistic 
7  is  the  pair  (r ,  y),  where  7  is  the  number  of  failures  observed  and  y  (defined 
by  Eq  5)  is  the  sum  of  the  erth  powers  of  the  operAtinft  times. 

The  next  subsection  briefly  discusses  a  particular  type  of  sampling  plan 
Experiment)  and  gives  the  distributions  of  7  and  7.  We  also  consider  the  time 
7  required  to  complete  the  experiment.  It  is  assumed  that  the  prior  distribu¬ 
tion  on  X  is  gamma-1  with  parameters  r  and  y\  Two  other  sampling  plans 
are  then  mentioned  in  the  subsection,  “Other  Sampling  Plans.* 


Test  n  Items  Simultaneously 
until  r  Failures  Occur 


without  Replacement, 


In  this  sampling  plan  n  items  are  placed  on  life  test  at  the  same  time  and 
the  experiment  is  terminated  when  r  of  them  have  failed.  In  the  special  case 
in  which  r  =  n,  the  failure  times  of  all  n  are  observed.  The  statistic  r  is  pre¬ 
determined  and  the  distribution  of  y  is  desired.  It  can  be  shown  (see  RAC-TP- 
2159)  that  this  distribution  is  inverted-beta-2: 


f(y\r',y'ir)  -  (lB.  <y| r.r'.yl  * 


y^1  (yT 


Note  that  this  distribution  does  not  depend  on  n.  The  time  7  to  complete  the 
experiment  (i.e.,  the  time  of  the  rth  failure)  does  depend  on  n;  the  density  func¬ 
tion  for  7  is  very  complicated,  but  the  pth  moment  is  given  by  (see  RAC-TP- 
2159): 

E<P|r\y';r.*)  -  E(P|l:r.ii)(yV/0r<r'-p/a)/m.  (10) 

Here  E (i*  |l;  r,  n)  is  the  expected  value  of  the  pth  power  of  the  rth  order  sta¬ 
tistic  in  a  sample  of  n  items.  Each  of  these  items  has  a  Weibull  lifetime  dis¬ 
tribution  with  X  =  1  and  the  same  value  of  or  used  in  the  Bayesian  analysis. 
The  first  moment,  E(?|l;r,  n),  is  tabulated  by  Harter10  for  a  ■  0.5(0. 5)4(l)8f 
r  =  1(1  Jn,  n  =  1(1)40.  The  pth  moment  may  be  obtained  by  replacing  or  by  a/p. 

Other  Sampling  Plans 

There  are  a  number  of  other  sampling  plans  that  may  be  considered  in 
practice  and  are  therefore  worthy  of  further  investigation.  Two  of  them  are 
as  follows: 

(1)  Test  items  sequentially  until  r  failures  occur.  As  soon  as  one  item 
fails,  another  one  is  placed  on  life  test.  The  statistic  r  is  fixed  and  y  and  7 
are  random  variables.  This  sampling  plan  is  discussed  in  RAC-TP-215.* 

(2)  Test  n  components ^simultaneously,  without  replacement,  and  stop  at 
time  t* »  min(Tw,  T),  where  isjhe  time  of  the  mth  failure  and  T  is  fixed  in 
advance.  In  this  case  7,  7,  and  ?  are  all  random  variables. 


5.  A  PROBLEM  IN  OPTIMA  r  SAMPLING 


To  illustrate  the  ,  ;ous  developments  we  present  an  analysis  (including 
prior,  posterior,  ami  pi  terior  analyses)  of  a  hypothetical  problem.  The 
sampling  plan  is  the  .  ;  the  one  considered  previously,  namely  testing  n 
items  simultaneously  v*  about  replacement,  until  r  failures  occur.  It  is  desired 
that  the  optimal  values  of  r  and  n  be  found. 

Problem  Statement 

A  manufacturer  must  supply  a  traveling- wave  tube  for  use  in  a  govern¬ 
ment  satellite.  He  can  either  use  one  of  type  T1  (act  a,)  whose  lifetime  dis¬ 
tribution  is  known  or  use  one  of  type  T2  (act  a2)  which  has  just  been  developed 
and  whose  lifetime  distribution  is  unknown.  He  is  willing  to  assume  that  the 
lifetime  distribution  of  type  T2  is  Weibull  with  a  particular  value  of  the  shape 
parameter  a.  He  assigns  a  gamma-1  prior  distribution  with  parameters  r  and 
V'  to  the  unknown  scale  parameter  A. 

The  manufacturer’s  terminal  utility  function  it,  (a.  A)  is  derived  as  follows. 

If  the  tube  used  in  the  satellite  functions  longer  than  time  ta  ,  benefits  with 
utility  k  accrue  to  the  manufacturer.  If  the  tube  fails  before  time  his  utility 
is  zero.  His  overall  utility  is  thus  k[Prob(tube  life  exceeds  t„  )],  and  therefore 

Ujiflj.A)  k  r~M°  kp  ul{a2.p), 

where  f 


If  A  is  a  random  variable,  so  is  p.  Since  the  lifetime  distribution  of  type  T1  Is 
known  we  have  «,(«,,  A)  =  uf  (a,,  p)  -  K,  a  constant. 

Prior  and  Posterior  Analysis 

It  is  convenient  in  this  problem  to  work  with  the  random  variable  p.  We 
assume  that  K  <  k  so  that  there  exists  a  breakeven  value  of  p,  called  pb ,  such 
that  u,(a, ,  pb)  =  u,  (a2,  pi)\  thus  p*  =  K/k.  Then  a,  is  preferred  for  p  <  ph,  and 
a,  is  preferred  for  p  >  ph . 

Since  p  is  a  random  variable,  if  terminal  action  *s  to  be  taken  on  the  basis 
of  prior  information  the  manufacturer  should  choose  an  act  o'  such  that 

CpK  fo'.pt  mmktrpU' ta.p)  . 

But  Ep  u,  (fl,,  p)  -  K  and  E^u,  (a,,  p)  ^  kp',  where  p  -  [y/(y  *■  tj)7'  is  the 
mean  of  the  prior  distribution  on  p.  Thus  the  optimal  act  under  the  prior  dis¬ 
tribution  is 

a.  fl,  ‘fe  -  A  • 

o,  <f  ?'  ’  ,v 

If  an  experiment  (life-test  sample)  on  tubes  of  type  T2  yields  an  outcome 
with  statistics  r  and  y,  the  posterior  distribution  on  A  will  be  gamma- 1  with 
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'  parameters  r"  =  t’  +  r  and  y"  =  y  +  y.  The  optimal  act  under  the  pos¬ 
terior  distribution  will  then  be 

*1  if  F  -  Pb ' 

where  p"  =  [y'/(y~  +  is  the  mean  of  the  posterior  distribution  on  p. 

Preposterior  Analysis 

In  considering  the  value  of  sampling  we  assume  that  sampling  and  termi¬ 
nal  utilities  are  additive,  so  that  u(e,  z  ,  a,  =  ut(a,  p)  -  cs(e,  z),  where 
cs  (e,  z)  is  the  cost  of  experiment  t  and  outcome  z.  Equation  8  then  becomes 

-  i if(e)  - 

We  have  E^jz  u,  (a,  p)  *  u,  (a,  p"),  where  a  tilde  is  placed  over  p~  to  indicate 
that  prior  to  the  experiment  it  is  a  random  variable.  We  next  note  that 

K  if  r**. 

y  ifr  > 

We  shall  consider  the  sampling  plan  treated  in  the  subsection  "Test  n 
Items  Simultaneously,  without  Replacement,  until  r  Failures  Occur.”  The 
statistic  7  is  a  random  variable  with  inverted-beta- 2  distribution  given  by  Eq 
9.  Since  p "  =  [y"/(y~  +  t®  )3r”,  we  obtain 

«;(<)  *  /^Kfi/3j(y|r.r'.yl<fy  *  JT  k  [<yyV<y+yM?>]r'f(/Jl  (y|r.r',y'My.  fl|) 
where  ^  is  defined  by 

Pi  *  [<yb*  >'V(yb  ♦  y’*  . 

Now  we  turn  tojhe  expected  cost  of  sampling  c£(r).  We  shall  assume  that 
c,(e,  7)  s  c,n  +  Cj?  ♦  Cjt*  ,  where  n  is  the  number  of  tubes  placed  on  life 
test,  7  is  the  duration  of  the  test,  and  the  ct  are  positive  constants.  We  then 
have 

cj(e)  -  c,n  ♦  CjE(iy.y';M«)  ♦  cs  Eff^r'.y'jr,*), 

where  E(t*  |r' ,  y'  ;  r ,  n)  is  given  by  Eq  10. 

The  final  step  in  preposterior  analysis  is  to  determine  the  optimal  sampling 
plan,  i.e.,  we  want  to  choose  r  and  n  to  maximize  u*(e)  =  u,#(f)  -  c*(e).  It  may 
happen,  of  course,  that  «*(«)  is  less  than  the  utility  of  immediate  terminal 
action  EpU,  (o',  p)  for  all  combinations  of  r  and  n;  in  this  case  no  experiment 
should  be  performed. 

Determination  of  the  Optimal  Sampling  Plan 

We  can  use  dynamic  programming  to  find  the  pair  (r*,  n#)  which  maxi¬ 
mizes  «*(«).  Let  n*(r,  *0  be  the  (expected)  utility  of  the  sampling  plan  t  =  (r,  »). 
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Define 


ti*(r)  max  u*(r.n)  . 

H 


Then 


u*<r*.rt*)  max  utr). 
r 

The  computations  required  to  obtain  W*(r)  are  relatively  minor  because  only 
c  *(c)  depends  upon  n;  u*(e)  depends  only  upon  r.  Determination  of  u*(r)  there¬ 
fore  amounts  to  minimizing  c*(r,  n)  with  respect  to  n  and  subtracting  the  re¬ 
sult  from  u*(r)  (we  write  u*(e)  =  uf(r)  when  c  specifies  r  failures).  Recall 
that  c*(r,n)  =  c,n  +  c2  E(t"| r'  y' ;  r,  n)  +  c,  E(7*  jr',  y’;  rf  n).  It  can  be  veri¬ 
fied  from  the  values  in  Harter10  that  E(tP  |r,  y' ;  r,  n)  is  a  convex  function 
of  n  for  fixed  r.  It  then  follows  that  c*(r,  n)  is  a  convex  function  of  n  for  fixed 
r  and  so  can  be  minimized  over  integer  values  of  n  fairly  readily. 

Calculation  of  u*(r)  from  Eq  11  involves  determination  of  values  of  the 
incomplete  beta-function.  Depending  on  the  computational  facilities  available, 
this  determination  might  be  done  only  for  selected  values  of  r.  Then  u#(r) 
could  be  computed  for  these  values  of  r,  and  the  results  graphed  as  a  function 
of  r.  The  graph  would  give  a  good  estimate  of  r*  and  this  estimate  could  be 
checked  by  computing  u*(r)  for  appropriate  values  of  r. 
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Appendix  A 


EXAMINATION  OF  A  PRIOR  DENSITY  FUNCTION 
WITH  TWO  PARAMETERS  UNKNOWN 


Here  we  shall  briefly  examine  a  prior  density  function  for  the  case  in 
which  the  parameters  A  and  q  of  the  Weibull  distribution  are  both  assumed  to 
be  unknown.  Although  this  density  function  possesses  property  4  of  Sec  2,  it 
does  not  really  have  property  5  and  definitely  does  not  possess  property  6. 

We  take  the  joint  density  function  on  (A,  a)  in  the  form  of  the  likelihood, 
expression  3: 

ju. sy. 4,. a'  /  <d,  J,  >  ’rvp  [,\  a  j.'J  Ky.rf,.  .J„>.  (12) 

where  0  s  A,  o  <  ■  r  >  1;  dt  >  0;  and 


’r  !'  ■  d,->  *  I  i  d;  »  '■  1  d. 


(13) 


Equation  12  defines  a  proper  density  function  if  some  d,  is  greater  than  1, 
since  this  suffices  to  show  that  k(r" ,  d, . dn)  is  finite. 

Now  if  the  prior  density  function  on  (A.  o)  is  f(A,  o\r  ,  d . dn  )  and 

the  likelihood  of  the  evidence  c  is  given  by  expi°ssion  3.  then  Bayes*  theorem 
yields 

((A.  j  f  •  r.d,.  ■*„’  (14) 

for  the  posterior  density  function.  Thus  the  prior  density  function  defined  by 
Eq  12  possesses  property  4.  Although  the  prior  and  posterior  density  functions 
have  similar  forms,  property  5  is  not  present  because  the  two  density  functions 
do  not  have  the  same  number  of  parameters.  Property  6  is  lacking  for  the  sim¬ 
ple  reason  that  the  constant  k(f  .  d,  ,  .  .  .  ,  dt|  )  cannot  be  evaluated  in  closed 
form. 
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